Astronomy & Astrophysics manuscript no. aal9999-12 


©ESO 2012 


October 8, 2012 





Starspot activity and rotation of tlie planet-liosting star Kepler-17 

A. S. Bonomo' and A. F. Lanza^ 



(N 

T— I 

O 

(N 
> 

o 

O; 
in 

pq 
^' 

o 



' INAF-Osservatorio Astrofisico di Torino, via Osservatorio, 20, 10025, Pino Torinese, Italy 

e-mail: bonomo@oato . inaf . it 
^ INAF-Osservatorio Astrofisico di Catania, Via S. Sofia, 78, 95123 Catania, Italy 

e-mail: nuccio . lanza@oact . inaf . it 



Received 12 July 2012 / Accepted 20 September 2012 



ABSTRACT 



Context. Kepler-17 is a G2V sun-like star accompanied by a transiting planet with a mass of « 2.5 Jupiter masses and an orbital 
period of 1.486 d, recently discovered by the Kepler space telescope. This star is highly interesting as a young solar analogue. It is 
also a good candidate for a test of the tidal theories for solar-like stars. 

Aims. We used about 500 days of high-precision, high-duty-cycle optical photometry collected by Kepler to study the rotation of the 
star and the evolution of its photospheric active regions. 

Methods. We applied a maximum-entropy light curve inversion technique to model the flux rotational modulation induced by active 
regions that consist of dark spots and bright solar-like faculae with a fixed area ratio. Their configuration was varied after a fixed time 
interval to take their evolution into account. Active regions were used as tracers to study stellar differential rotation, and planetary 
occultations were used to constrain the latitude of some spots. 

Results. Our modelling approach reproduces the light variations of Kepler-17 with a standard deviation of the residuals comparable 
with the precision of Kepler photometry. We find several active longitudes where individual active regions appear, evolve, and decay 
with lifetimes comparable to those observed in the Sun, although the star has a spotted area « 10-15 times larger than the Sun at 
the maximum of the 11 -yr cycle. Kepler-17 shows a solar-like latitudinal differential rotation, but the fast spot evolution prevents a 
precise determination of its amplitude. Moveover, the stars shows a cyclic variation of the starspot area with a period of 47.1 ± 4.5 d, 
particularly evident during the last 200 days of the observations, similar to the solar Rieger cycles. Possible effects of the close-in 
massive planet on stellar photospheric activity cannot be excluded, but require a long-term monitoring to be unambiguously detected. 
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1. Introduction 

The Kepler space mission has recently announced the discov- 
ery of the hot Jupiter Kepl er- 17b. It orbits a G2V star with 
= 14.14 in P = 1.486 d (iDesertet al.ll201lL hereafter DI 1). 
Thanks to additional spectra acquired with the SOPHIE spec- 
trograp h at the Observatoire de Haute Provence. iBonomo et al. I 
(1201 2h (hereafter B 12) slightly refined the characterization of the 
Kepler-17 planetary system. In particular, they found the host 
star to be slightly hotter and younger than D 1 1 , with an effective 
temperature of 578 1 + 85 K and an age < 1 .8 Gyr. 

The hosting star Kepler-17 is remarkably active, as shown by 
the out-of-transit variations in the Kepler light curve with a peak- 
to-peak amplitude of * 4%. They are produced by the rotational 
modulation of active regions, starspots and faculae, on the stellar 
disc. By means of a periodogram analysis of the light curve, DI 1 
estimated the stellar rotation period to be Prot - 1 1.89 ± 0.15 d. 
Intriguingly, this value of Prot is an integer multiple of the orbital 
period, specifically eight times. This could be just a coincidence 
or the result of a star-planet interaction, as pointed out by D 1 1 . 

In addition, the planet occults starspots during transits, giv- 
ing rise to the typ ical bumps observed in the bottom of the tran- 
sit pr ofiles (e.g., ISilva-Valio et al. 1 l20lOHSilva -Valio & Lanza I 
12011 1). The integer ratio between P^ot and P allowed DI 1 to see 
a "stroboscopic" effect with the short-cadence Kepler data: the 



Kepler magnitude (cf. Sect.[2]l. 



spots are "mapped" by the planet each 45 ° in longitude. The 
study of these spot-crossing events reveals that the planet's or- 
bit is prograde and the projected spin-orbit angle is smaller than 
10° - 15° (see Dll for more details). 

Kepler-17 is therefore an excellent candidate to study the 
magnetic activity in a solar-like star younger than the Sun, 
through the modelling of its out-of-transit flux variations. A 
similar analysis has been performed for several CoRoT planet- 
hosting stars. One of them is CoRoT-2, an active G7V star with 
a rotation period o f 4.5 d and an age younger than 0.5 Gyr 
(lAlonso et al.ll200i ). The spot-modelling of its light curve re- 
vealed two active longitudes and a cyclic oscillation of the to- 
tal spotted area wit h a period of 28.3 + 4.3 d, similarly to the 
solar Rieger cycles dLanza et al ]L2009a'). Investigating the mag- 
netic activity of G dwarfs with different ages woul d in princi- 
ple allow one to study the Sun i n time. Recently, Frasc a et al. I 
( 2011) and lFrohlich et al.l (1201 2h have applied a spot-modelling 
to Kepler photometric time series to study three young G- and K- 
type dwarf stars whose ages range between » 50 and ^ 200 Myr, 
deriving information on their spot evolution and surface differ- 
ential rotation. 

Owing to the faintness and low projected rotation velocity of 
Kepler-17 (vsin; ^ 5 km s"'), an investigation of its magnetic 
activity through Doppler imaging is not feasible. Therefore, 
modelling the stellar variability and studying the distortions of 
the transit profile when the planet occults starspots are the only 
techniques that allow us to derive active longitudes where spots 
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preferentially form, determine the lifetime of active regions and 
activity complexes, estimate a minimum amplitude for the stel- 
lar differential rotation, and possibly discover short-term and/or 
long-term activity cycles, depending on the length of the time 
series. 

2. Observations 

The Kepler space telescope has an aperture of 95 cm and is de- 
signed to yield nearly continuous, high-precision photometry in 
the passband 423 - 897 nm for a: 150, 000 stars i n a fixed field 
of view to search for pl anetary transits (see, e.g., iBorucki et al.l 
l2010t iKoch et all 12010) . It orbits the Sun on an Earth-traihng 
orbit and, to keep the solar arrays pointed towards the Sun, a 
re-orientation of the spacecraft is required every « 90 days, a 
time interval called "a quarter" in the Kepler jargon. The re- 
orientation of the telescope produces an offset in the photon 
counting for a given star because its photometric mask is re- 
defined on a different CCD in the focal plane. 

The light curve of Kepler- 17, publicly available at the MAST 
archiveQ, covers more than fifteen months of photometric mea- 
surements, from 2009 May 13 to 2010 August 23. Observations 
are distributed along the six quarters Q1-Q6. The raw data 
with the long-cadence temporal sampling, i.e. one point each 
29.42 min, was used for our work. Short-cadence data (one point 
per minute) available for the last three quarters (Q4-Q6) are not 
particularly useful for our purpose because we are interested in 
modelling the out-of-transit variations on the timescale of stellar 
rotation or longer. 

Because of the flux offsets between adjacent quarters and 
the long-term instr umental trends within each quarter (see 
iJenkins et al.ll2010al and Fig. 1 in Dl 1 for our specific case), the 
light curves corresponding to diff'erent quarters were separately 
treated as follows. First, planetary transits were removed from 
each quarter The flux contamination due to starfield crowding, 
as estimated by the Kepler team, was subtracted from the me- 
dian value of the flu x. Steep variations after the safe modes 
(iJenkins et al. Il2010b^ were removed and long-term trends of 
clear instrumental origin were coiTected by fitting a parabola. 
Lastly, the flux of each quarter was normalized to its median 
value then nearly matched the endpoints of adjacent quarters. 
The final light curve obtained by combining the six quarters con- 
tains 20924 data points and is shown in Sect. 15. II (see Fig. |2]i. 
The median of the eiTors of the single photometric measure- 
ments is 2.18 X 10""* in relative flux units. 



3. Spot modelling of wide-band light curves 

Reconstructing the surface brightness distribution from the ro- 
tational modulation of the stellar flux is an ill-posed problem, 
because the variation of the flux vs. rotational phase contains 
information only on the distribution of the brightness inhomo- 
geneities vs. longitude. The integration over the stellar disc ef- 
fectively cancels any latitudinal information, particularly when 
the inclination of the rotation axis along the line of sight is 
close t o 90°, as in the present case (see Sect. |4]and iLanza et all 
I2009ah . Therefore, we need to include a priori information in 
the Ught curve inversion process to obtain a unique and stable 
map. This is done by computing a maximum entropy (here- 
inafter ME) map, which has been proven to best reproduce 
active region dist ribution and area variations for the Sun (cf. 
ILanza et al.ll2007l) . For a different modelling approach, based on 



discrete starspots, see, e.g.. lMosser et alJ (l2009l) . A comparison 
with other modelling approaches is given in, e.g.. Frohlic h et al] 
(2009); Huberet al. (2010); Silva-Valio & Lanza (2011). 

In the present model, the stellar surface is subdivided into 
elements, i.e., into 200 squares of side 18°, with each element 
containing unperturbed photosphere, dark spots, and facular ar- 
eas. The fraction of the fc-th element covered by dark spots is 
indicated by its filling factor fk, the fractional area of the faculae 
is Qfk, and the fractional area of the unperturbed photosphere is 
I - {Q + l)fk- The contribution to the stellar flux coming from 
the k-ih surface element at time tj, where j - 1 , is an index 
numbering the points along the light curve, is given by 



^Fkj = h{txkj){i - (2 + m + cjk+ 

+cf(i - ^lkj)]\Ak^lkjw{mi)^ 



(1) 



where /q is the specific intensity in the continuum of the unper- 
turbed photosphere at the isophotal wavelength of the observa- 
tions, Cs and Cf are the spot and facular contrasts, respectively (cf. 
ILanza et a n i2004 . A,, is the area of the A:-th surface element. 



MMkj) 



1 ifpkj>0 
if Pkj < 



(2) 



is its visibility, and 



pitj = cos iffkj = sin / sin Ok cos[4 + Q.(tj - to)] + cos / cos Ok, (3) 

is the cosine of the angle i/'^j between the normal to the sur- 
face element and the direction of the observer, with / being the 
inclination of the stellar rotation axis to the line of sight, 0^ 
the colatitude and 4 the longitude of the ^-th surface element, 
Q. the angular velocity of rotation of the star (Q = In/Proi), 
and to the initial time. The specific intensity in the continuum 
varies accordin g to a quadratic limb-darkening law, as adopted 
by Lanza et al. (2003j for the Sun, viz. /q cc ap + bpp + Cpfi^. The 
stellar flux computed at time fj is then: F(tj) = Yjk ^Fkj- To war- 
rant a relative precision of about 10"^ in the computation of the 
flux F, each surface element is further subdivided into elements 
of 1° X 1° and their contributions, calculated according to Eq. ([T]i, 
are summed up at each given time to compute the contribution 
of the 18° X 18° surface element to which they belong. 

We fitted the light curve by changing the value of the spot- 
filling factor / over the surface of the star while Q was held 
constant. Even fixing the rotation period, the inclina tion, and 
the spot and facular contrasts (see iLanza et"aL I I2OO7I for de- 
tails), the model has 200 free parameters and suffers from non- 
uniqueness and instability. To find a unique and stable spot 
map, we applied ME regularization, as described in Lanza et alj 
(200_2), by minimizing a functional Z, which is a linear combi- 
nation of tile;^f^ and the entropy functional S ; i.e.. 



(4) 



where / is the vector of the filling factors of the surface ele- 
ments, A > a Lagrangian multiplier determining the trade-off 
between light curve fitting and regularization, and the expres- 
sion of S is 



Wk 



.Alog- + (l 
m 



•/t)log 



l-/k 



1 



(5) 



http://archive.stsci.edu/k;epler/data_search/search.php 



where is the relative area of the A:-th surface element (total 
surface area of the star =1) and m the default spot-covering fac- 
tor that fixes the limiting values for ft as: m < < (1 - m) . 
In our modelling we adopted m - 10"^ (cf. iLanza et"anil998h . 
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The entropy functional S is constructed in such a way that it at- 
tains its maximum value when the star is virtually immaculate 
(fk - m for every k). Therefore, by increasing the Lagrangian 
multiplier A, the weight of S in the model is increased and the 
area of the spots is progressively reduced. This gives rise to sys- 
tematically negative residuals between the observations and the 
best-fit model when A > 0. 

The optimal value of A depends on the information content 
of the light curve, which in turn depends on the ratio of the am- 
plitude of its rotational modulation to the average standard de- 
viation of its points. In the case of Kepler-17, the amplitude of 
the rotational modulation is x 0.044, while the nominal standard 
deviation of the points is s; 2.18 x 10""* in relative flux units (see 
Sect.|2]i, giving a signal-to-noise ratio of * 200. 

To fix the optimal value of the Lagrangian multiplier A, we 
compared the modulus of the mean of the residuals of the regu- 
larized best fit ljUregl with the standard error of the residuals them- 
selves, i.e., eo = (To/ VA^, where ctq is the standard deviation 
of the residuals of the unregularized best fit and A^ the number 
of points in each fitted subset of the light curve having a dura- 
tion Aff (see below). We iterated until l/iregl - fo (cf. the case of 
CoRoT-2 in iLanza et alj|2009al) . 

For the Sun, by assuming a fixed distribution of the filling 
factor, it is possible to obtain a good fit of the irradiance changes 
only for a limited time interval Aff, not exceeding 14 days, which 
is the lifetime of the largest sunspot groups dominating the irra- 
diance variation dLanza et al.ll2003l ). For other active stars, the 
value of Aff must be determined from the observations them- 
selves, looking for the maximum data extension that allows us a 
good fit with the applied model (see Sect.lU. 

The optimal values of the spot and facular contrasts and of 
the facular-to-spotted area ratio Q in stellar active regions are 
unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product CfQ, so we can fix Cf and vary 
Q, estimating its best value by minimizing the of the model, 
as shown in Sect. |4] Since the number of free parameters of the 
ME model is large, we used the model of Lanza et al. (2003) 
to fix the value of Q. It fits the light curve by assuming only 
three active regions to model the rotational modulation of the 
flux plus a uniformly distributed background to account for the 
variations of the mean light level. This procedure is the same as 
ado pted for, e.g., CoRoT-2 and CoRoT-4 to fix the value of Q 
(cf. iLanza et al]|2009al; iLanza etaL 1 l2009bl) . 

We assumed an inclination of the rotation axis of Kepler-17 
of / = 87.2° (see Sect. |4]i. Since the information on spot lati- 
tudes that can be extracted from the rotational modulation of the 
flux for such a high inclination is negligible, the ME regular- 
ization virtually puts all spots at the sub-observer latitude (i.e., 
90° - / ~ 3°) to minimize their area and maximize the entropy. 
Therefore, we are limited to mapping only the distribution of the 
active regions vs. longitude, whi ch can be achieved with a reso- 
lution of at least ^ 40° - 50° (cf. iLanza et al.ll2007l[2009 a). Our 
ignorance of the true facular contribution to the light modulation 
may lead to systematic errors in the active region longitudes de- 
rived by our model, as discussed by Lanza et al.. (.2007.) for the 
Sun. 



4. Model parameters 

The basic stellar parameters are taken from D 11 and B 1 2 and are 
listed in Table [T] The limb-darkening parameters in the Kepler 
bandpass were derived from the model of the transit as discussed 
inDll. 



Table 1. Parameters adopted for the modelling of the light curve 
of Kepler-17. 



Parameter Value Ref." 



Star mass (Mq) 


1.16 


B12 


Star radius (Rq) 


1.05 


B12 


Teff (K) 


5780 


B12 


log g (cm s-2) 


4.53 


B12 




0.333 


BL12 


b. 


0.929 


BL12 




-0.262 


BL12 


-Plot (days) 


12.01 


BL12 


^rot 


4.66 X 10-^ 


BL12 


Inclination (deg) 


87.22 


Dll 




0.677 


L04 


Cf 


0.115 


L04 


Q 


1.6 


BL12 


Aff (days) 


8.733 


BL12 



" References: B 12:lBonomo et al. I (f2012): Dl l: lDesert et all ( 1201 ih : 
L04: lLanzaetalJ ( l2004i) : BL12: present study. 



The rotation period for our spot modelling was initially fixed 
at exactly eight orbital periods of the planet, following Dll, who 
found that there were spots occulted by the planet that reap- 
peared at the same phase after that time interval, i.e., 11.89 d. 
However, we found that the migration rates of the other spots 
revealed by the modelling of the out-of-transit photometry was 
minimized by assuming a slightly different rotation period, i.e., 
12.01 d, which we adopted for our modelling. 

The polar flattening of the star owing to the centrifugal po- 
tential was computed in the Roche approximation with a rotation 
period of 12.01 d. The relative difference between the equatorial 
and the polar radii is erot - 4.66 x 10"^, which induces a com- 
pletely negligible relative flux variation of » 2 x 10"^ for a spot 
coverage of » 4 percent, as a consequence of the gravity darken- 
ing of the equatorial regions of the star. 

The inclination of the stellar rotation axis is constrained by 
the observation of Dl 1 that the light anomalies detected in suc- 
cessive transits are compatible with the same spots being repeat- 
edly occulted as they move over the stellar disc owing to stellar 
rotation. This is possible only if the disc-projected trajectory of 
the spots remains inside the belt occulted by the planet. In other 
words, this is an indication that the sky-projected misalignment 
between the stellar spin and the orbital angular momentum is 
smaller than ±(10° - 15°) (cf. Dll for details). Although the 
possibility that the stellar spin is not aligned with the orbital 
angular momentum cannot be excluded with certainty because 
only the sky-projected angle between the two vectors has been 
measured, we regard this occurrence as highly improbable and 
assume that the rotation axis is inclined with respect to the line 
of sight of the same angle as the orbit normal. Moreover, the 
stellar v sin / = 4.7 + 1 .0 km s"' and the estimated stellar radius 
provide an inclination close to 90° for the stellar rotation period 
(cf. Dll). 

The maximum time interval Aff that our model can accu- 
rately fit with a fixed distribution of active regions was deter- 
mined by dividing the total interval, T = 497.785 d, into A^f 
equal segments, i.e., Aff = T/Nf, and looking for the minimum 
value of A^f that allows us a good fit of the light curve, as mea- 
sured by the;t'^ statistics. We found that for A^f < 57 the quality 
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of the best fit degrades appreciably with respect to higher values, 
owing to a substantial evolution of the pattern of surface bright- 
ness inhomogeneities. Therefore, we adopted A^f = 57, so that 
Aff - 8.733 d is the maximum time interval to be fitted with a 
fixed distribution of surface active regions to estimate the best 
value of the parameter Q (see below). 

To evaluate the spot contrast, we adopted the same mean 
temperature difference as deri ved for sunspot g roups from their 
bolometric contrast, i.e. 560 K dChapman et alll 994). The effec- 
tive temperature of the unspotted photosphere is 5780±85 K, i.e., 
very similar to that of the Sun (cf. B12). In other words, we as- 
sumed a spot effective temperature of 5220 K, yielding a cont rast 
Cs - 0.677 in the bolometric passband (cf. iLanza et"al]|2007h . A 
different spot contrast changes the absolute spot coverages, but 
neither significantly affects their longitudes or their time evolu- 
tion, as discussed in detail by Lanza et al. (2009a). Therefore, 
since the spot temperature is only estimated, we neglected the 
difference in the contrast between the Kepler bandpass and 
the bolometric passband in our analysis. 

In our model, the facular contrast is assumed to be solar-like 
with Cf - 0.115 ( Lanza et al. 2004) . The best value of the area 
ratio Q between the faculae and the spots in the a ctive regions 
was es timated by means of the three-spot model bv lLanza et alj 
( I2003L cf. Sect. [2). In Fig. □ we plot the ratio X^lxlai^ of the 
total of the composite best fit of the entire time series to its 
minimum value versus Q, and indicate the 9 5 percent confi- 
denc e level as derived from the F-statistics (e.g., iLampton et alj 
[1976). The choice of Aff = 8.733 d allows us to fit the rotational 
modulation of the active regions for the longest time interval dur- 
ing which they remain stable, modelling both the flux increase 
due to the facular component when an active region is close to 
the limb and the flux decrease due to the dark spots when the 
same region transits across the central meridian of the disc. In 
this way, a measure of the relative facular and spot contributions 
can be obtained, leading to an estimate of Q. Fig. [1] shows that 
the best value is, Q - 1.6, with an acceptable range extending 
from X 0.6 to a: 2.6. Therefore, we adopted Q - 1.6 for our 
modelling in Sect. |5] We comment on the value of Q in more 
detail in Sect.|6] 




12 3 4 

Q 



Fig. 1. Ratio of the of the composite best fit of the entire 
time series of Kepler-17 to its minimum value vs. the parameter 
Q, i.e., the ratio of the area of the faculae to that of the cool 
spots in active regions. The horizontal dashed line indicates the 
95 percent confidence level for;f^/A'^in' determining the interval 
of acceptable Q values. 



5. Results 

5. 1 . Light curve models 

We applied the model of Sect. [3] to the out-of-transit light curve 
of Kepler-17, considering time intervals Aff - 8.733 d. The best 
fit without regularization {A - 0) has a mean yUi-es = 3.845 x 10"^ 
and a standard deviation of the residuals ctq = 3.032 x 10""* in 
relative flux units. The Lagrangian multiplier A is iteratively ad- 
justed until the mean of the residuals jUres - -1-583 x 10"^ - 
(To/ -s/N, where A^ - 367 is the mean number of points in each 
fitted light curve interval Aff; the standard deviation of the resid- 
uals of the regularized best fit is cr = 3.31 1 x 10 

The composite best fit to the entire light curve is shown in the 
upper panel of Fig.|2] while the residuals are plotted in the lower 
panel. The best fit is always very good, with a standard devia- 
tion of the residuals x 1 .52 times the median of the eiTors of the 
photometric measurements as given by the Kepler pipeline. The 
distribution of the residuals is plotted in Fig.[3]and is well fitted 
by a Gaussian with a standard deviation of 3.056 x 10"'* for abso- 
lute values of the residuals lower than ^ 6x 10 in relative flux 
units. For residuals greater than a; 6 x 10 '* there is a remarkable 
asymmetry in the distribution with an excess of positive resid- 
uals. A periodogram of the residuals shows a highly significant 
peak at a period of 1 .505 ± 0.005 d, very close to the orbital pe- 
riod, and several lower peaks at periods between a 1 .2 and a 3 d. 
Putting the residuals in phase with the orbital period, we see the 
occultation of the planet and a possible phase-dependence of the 
reflected light (see B12). 

5.2. Longitude distribution of active regions and stellar 
differential rotation 

The distribution of the spot-filling factor / vs. the longitude and 
the time is plotted in Fig. |4] The longitude zero corresponds to 
the point intercepted on the stellar photosphere by the line of 
sight to the centre of the star at BJD 2454964.5 109, i.e., the sub- 
observer point at the initial epoch. The reference frame rotates 
with the star with a fixed period of 12.01 d and the longitude 
increases in the same direction as the stellar rotation and the or- 
bital motion of the planet. This is consisten t with the reference 
frames adopted in our previous studies (e.g., iLanza et alJ2009al : 




e;= 0.95 



-0.002 -0.001 0.000 0.001 0.002 

Residuals 

Fig. 3. Distribution of the residuals to the composite ME- 
regularized best fit shown in Fig. [2] The dashed line is a Gaussian 
fit to the residual distribution, while the vertical dotted line 
marks the zero value of the residuals. 
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100 



200 300 

Time (BJD - 2454900) 



400 
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Fig. 2. Upper panel: The ME-regularized composite best fit to the out-of-transit light curve of Kepler-17 obtained for Q - 1 .6. The 
flux is in relative units, i.e., normalized to the maximum observed flux along the light curve. Lower panel: The residuals from the 
composite best fit versus the time. 



iLanza et"al 1 l2009bh . but does not allow a direct comparison of 
the mapped active regions with the dips in the light curve. 

Our map shows that the individual starspots evolve with 
timescales of tens of days, which makes it difficult to trace the 
evolution and migration of the active regions in an unambiguous 
way. Nevertheless, two main active longitudes, where individ- 
ual spots form and evolve, can be identified with confidence and 
appear to rotate on the whole with the rotation period of our ref- 
erence frame. The migration rate of individual spots or groups of 
spots within each longitude is variable, which suggests that indi- 
vidual spots are forming at different latitudes on a differentially 
rotating star or, alternatively, there are several spots at close lon- 
gitudes that are evolving to mimic spot migration. Unfortunately, 
there is no information on the starspot latitude from our map- 
ping technique and even the hemisphere cannot be determined 
because the inclination is close to 90° (cf. Sect.O. Nevertheless, 
some additional information can be extracted in this particular 
case from the occultation of the spots by the planet during tran- 
sits. Specifically, we can exploit the results of Dll to constrain 
the latitude of some of the starspot trails that are seen in Fig.lD 

iDesert et alJ (|201 ll) identified five starspots that are repeat- 
edly occulted every eight transits and labelled them A, B, C, D, 
and E (cf. their Figs. 11 and 12). The transit profile distortions 
associated with spot C are barely visible in their Fig. 11, thus 
spot C is certainly significantly smaller than the other four spots. 
We plot in Fig. |4] the migration of the starspots detected by Dl 1. 
The initial epoch £(0) of Dll is equal to the epoch of the first 



mid transit as reported in their Table 3 (Desert, priv. comm.) so 
we can easily convert their phases of maximum starspot visibil- 
ity into longitudes in our reference frame. We find a good asso- 
ciation between their starspots A, B, D, and E and some trails of 
spots found with our approach. For starspot C we find almost no 
coincidence, as expected owing to its smaller spotted area. The 
resolution of the t ransit-mapp ing met hod is of the order of some 
degre es (cf., e.g., ISilva-Valio et al. 20Tol;ISilva-Valio & Lanzal 
1201 iL and references therein) while that of our mapping based 
on the out-of-transit light curve is approximately 40° - 50°. 
Therefore, we are not able to resolve the occulted spots with 
the same detail as Dll, especially when there are several small 
spots close in longitude. Note that a detailed comparison of the 
spot maps based on the in-transit and out-of-transit photometry 
is outside the scope of the present work because it requires a 
detailed modelling of th e transit profile distortions i nduced by 
occulted starspots (see iSilva-Vaho & Lanza 1 1201 iL for such 
a detailed comparison for CoRoT-2). Here we limit ourselves 
to exploiting the preliminary modelling of D 1 1 to find that our 
trails of spots with a positive migration rate, i. e., rotating faster 
than our reference frame, are located inside the belt that is oc- 
culted by the planet. A faster migration is particularly evident in 
the case of our spot trails associated with their spots B and D, 
at least for a significant part of the considered time interval. On 
the other hand, the intrinsic evolution of the pattern makes the 
determination of the migration rates of our spots associated with 
their spots A and E less certain, although a positive migration 
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Fig. 4. Distribution of the spot-filling factor vs. longitude and time for Q - 1.6. The values of the filling factor were normalized 
to their maximum /Inax - 0.01553 with orange-yellow indicating the maximum and dark blue the minimum (see the colour scale 
on the right lower corner of the figure for the correspondence between the colour and the normalized filling factor). Note that the 
longitude scale is extended beyond 0° and 360° to help following the migration of the starspots. The tracks of the five spots occulted 
during the planetary transits after Dll are also reported. The open circles mark the time intervals of eight transits after which the 
same spots are detected again during the transits. The straight lines connecting the circles trace the migration of those spots in our 
reference frame; each line is labelled with the name of the corresponding spot, as indicated in Fig. 11 of Dl 1. 



is suggested. On the whole, these results indicate that the oc- 
culted spots located at low latitudes are rotating faster than the 
overall pattern of the active longitudes, thus giving evidence of a 
solar-like differential rotation in Kepler-17, i.e., with the equator 
rotating faster than the poles. 

Finally, we note that the time resolution of our spot models 
is not adequate to look for a possible modulation of the stellar 
activity with the orbital period of i ts c lose-in r nassiv e planet, as 
suggested bv lShkolnik et all (l2008h or lLanzal (1201 ll) . Given the 
short orbital period of the planet, a different approach should 
be used to search for signatures of a possible star-planet in- 
teraction using sho rt-cadence data, as for CoRoT-2 (cf., e.g., 
iPagano et al.1l2009t) . 

5.3. Variation of tlie spotted area 

The lifetimes of the active longitudes traced in Fig.|4]range from 
ss 75 to X! 330 days and possibly longer, given the limited time 
interval covered by our observations. Individual active regions 
have a lifetime of about a few tens of days. This is similar to what 
we observe in the Sun where complexes of activity consisting of 
several active regions, forming approximately around the same 
longitude, have lifetimes up to 5 - 6 months. The duration of the 



light anomalies observed during the transits reaches to » 0.02 
in phase units, indicating a size of the largest occulted active 
regions of at least x 40° in longitude, i.e., 3-4 times that of the 
largest sunspot groups. 

The variation of the total starspot area vs. the time is plotted 
in Fig. |5] The error bars have a semi-amplitude of three stan- 
dard deviations as derived from the uncertainty of the photome- 
try, but systematic errors associated with the assumptions of the 
data reduction and the model are not included. Specifically, the 
variations on time scales longer than the duration of one quarter, 
i.e., 90 - 100 days, cannot be reconstructed from Kepler pho- 
tometry because of the flux jumps from one quarter to the next 
(see Sect. |2]). Therefore, we are limited to study the variations 
on timescales shorter than three months. Moreover, the pres- 
ence of gaps in the observations can introduce systematic errors 
in the measurement of the spotted area. This is the case of the 
gap beginning at BJD 2454997.982 for a duration of 5.0268 d. 
Since the time interval adopted for our individual best fits is 
Aff = 8.733 d, this loss of data implies a systematically lower 
value of the total spotted area because the ME regularization re- 
moves spots at the longitudes not constrained by the observa- 
tions. Fortunately, the other gaps in the photometric time series 
are much shorter than Aff, thus no other value of the area is sig- 
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nificantly affected. To show the distribution of the gaps, we plot 
line segments with a length equal to their duration at the level 
0.04 in Fig.|5]considering all interruptions with a duration longer 
than 24 hours. In the analysis presented below, we discarded the 
area value at BJD 2455003.81 because it is affected by the gap 
that started at 245997.982. 

The Lomb-Scargle periodogram of the entire time series of 
the area values is plotted in Fig. |6] (solid line) together with the 
power level correspondin g to a false-alarm p robability of 0.01 
as evaluated according to iHorne & Baliunas ] ll986). The main 
peak corresponds to a period of 47.1 ±4.5 d and has a false-alarm 
probability (PAP) of 2.0 percent. The value of the PAP was con- 
firmed by performing an analysis of 50,000 random Gaussian 
noise time series with the same sampling as our area data se- 
ries. We also plot the periodogram of the time interval from 
BJD 2455230.869 to 2455457.929 (dashed Une) because we see 
a regular oscillation of the total spotted area with a period of 
Si 48 d during that interval in Pig. |5] The main periodogram 
peak corresponds to a period of 48.2 ± 9.0 d with an PAP of 0.44 
percent. 

The time variation of the frequency of the spotted area mod- 
ulation is best represented by means of a wavelet amplitude. We 
plot in Pig. Q the amplitude of the Morlet wavelet versus the 
period and the time (see, e.g., iHempelmann & Donahue 1 [T997[ 
for details). The wavelet parameters are adjusted to give a time 
resolution of » 100 d for a period of about 48 d and a relative 
period resolution of APjP x 0.06. We see that in the initial part 
of the dataset there is a periodicity of a; 30 d that corresponds 
to secondary peaks in the periodogram of the whole time se- 
ries whose false-alarm probability is > 50 percent (cf. Pig. |6] 
where the frequency resolution is better than in the case of the 
wavelet). On the other hand, during the second half of the time 
series we see a clear periodicity of » 50 d, which corresponds 
to the significant peak in the periodogram. We conclude that the 
total spotted area of Kepler- 17 showed an oscillation with a pe- 
riod of 47.1 + 4.5 d. This behaviour is reminiscent of the short- 
term oscillations of the total sunspot area found close to the 
maxima of some of the 11 -year solar cycles. They were called 
Rieger cycles because they were first dete cted in the perio dic - 
ity of occurrence of large solar flares by Rieger et al. I (Il984l) . In 
the Sun, the periodicity is ^ 160 d with small variations from 
one cycle to the other, although only some of the sunspot max- 
ima show evidence of this short-term periodicity (Oliver et al. ' 
[1998; Krivova & Solanki 2002; Zaqarashvili et al. 2010). A be- 
haviour similar to that of Kepler- 17 was found in CoRoT-2, a 
G7V star that showed a Rieger- type cycle in the variation of its 
spotted area with a period of 28.9 + 4.3 d dLanza et al.ll2009al) . 

6. Discussion and conclusions 

The application of a spot model to reproduce the optical light 
curve of Kepler- 17 shows that the spot pattern is almost stable 
for a timescale of a: 9 d because the residuals to our best fits 
have a standard deviation of cr ~ 3.3 x 10 in relative flux units 
that is only » 50 percent greater than the mean error attributed 
to the photometric measurements by the Kepler pipeline. Our 
facular-to-spotted area ratio Q - 1 .6 is significantly lower than 
the value Qq = 9.0 ado pted for the modelh ng of the light curves 
of the Sun-as-a-star bv iLanza et al.l (l2007h . However, this lower 
value of Q i s typical of sun-like stars that are more activ e than 
the Sun (cf. Lanza et al.°200 9aLlLanza et al. 1120 1 01120 lla llbS and 
suggests an increasing weight of the dark spots in the photomet- 
ric variations a s a star becomes mo re a ctive, as indicated by th e 
results of, e.g. jRadick et al.l(ll998h and lLockwood et al.l(l2007h . 




00 
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Time (BJD-2454900) 

Fig. 5. Total spotted area as derived from the regularized ME 
models vs. time for Q = 1.6. The lower horizontal ticks mark 
the gaps in the photometric time series longer than 24 hours. 



FAP = 0.01 




Period (d) 

Fig. 6. Lomb-Scargle periodogram of the variation of the spotted 
area. The solid line gives the normalized power vs. the period for 
the whole time interval, while the dashed line gives the power for 
the time interval from BJD 2455230.8696 to BJD 2455457.9294 
with the same normalization as adopted for the periodogram of 
the whole interval. The horizontal dotted line marks the 99 per- 
cent confidence level (PAP = 0.01). 

Our models found several active longitudes that clustered on 
opposite hemispheres with a separation of a 180° for more than 
half the observation interval. This explains the two peaks found 
in the periodogram of the light curve by Dl 1, one corresponding 
to t he rotation peri o d and t he other to its first harmonic. 

iBonomo et al. ] (120 12h suggested an age younger than 
1.8 Gyr, while Dl 1 derived an age of 3.0 ± 1.6 Gyr for Kepler- 
17. As noted by B12, the age determined by means of standard 
gyrochronology is only 0.9 + 0.2 Gyr, while considering the ef- 
fects of the close-in plane t on the evolution of the stellar angular 
momentum ( Lanza l l20T^. an age of 1.7 ± 0.3 Gyr is estimated 
that seems more compatible with the age found from isochrone 
fitting. 

Since the rotation period of the star is longer than the orbital 
period, tides remove angular momentum from the orbit to spin 
up the star and lead to orbital decay. The timescale for the engulf- 
ment of the planet can be estimated according to Ogilvie & Lin] 
(120071) as Ta ^ 0.048(es/10*) Gyr, where is the modified 
tidal quality factor of the star This timescale is much shorter 
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than the lifetime of the system on the main sequence if we adopt 
2s ^ 10^, i.e., the value derived from the observed circulariza- 
tion periods of close binary systems in open clusters of differ- 
ent ages. Together with the observations of several other stars 
with massive planets on very tight orbits, this suggests that 
should be much higher (i.e., the tidal dissipation much lower) in 
those star-planet systems than in close binary systems that con- 
sist of two main-sequence stars. The difference in the value 
could arise because the stellar rotation is far from being syn- 
chronized with the orbital motion of its planet. In this case, con- 
sidering the dissipa tion of the tides inside the convection zone, 
lOgilvie & Liiil (120071) predicted > 5 x 10'', which impHes an 
infall timescale longer than the system lifetime. 

A lower limit on the value of can be set by an accurate 
timing of the transits over a time interval of a few decades be- 
cause for 2s = 10^ we expect a variation of the orbital period 
of APorb/^'orb ~ 5 X 10"** in ten years. It produces a variation 
of a; 8 s in the epoch of the mid transit in 10 years. The ac- 
curacy reported by D 11 is + 2.4 s for their initial transit epoch, 
implying that a of about 10^ should give an orbital period 
acceleration detectable in a few years with a space-borne pho- 
tometer. The recently approved extension of the Kepler mission 
till 2016 is therefore an interesting opportunity to perform such 
measurements. A model of the light perturbations that are caused 
by the spots occulted during the transits may possibly be used to 
improve the accuracy since Dll found small O - C timing os- 
cillations with a period of approximately half the stellar rotation 
period, i.e., likely associated with starspots on opposite stellar 
hemispheres (cf. their Fig. 10). 

We can estimate an approximate lower limit to the differen- 
tial rotation of Kepler-17 finding AQ/O » 0.10 - 0.16 from the 
migration rates of the different trails of spots as seen in Fig. |4] 
Given the rapid evolution of the individual spots, this value is 
not only approximate, but depends critically on the way indi- 
vidual spot trails are identified. Therefore, we caution to regard 
this estimate as remarkably uncertain (cf. Sect. I5.2l i. For com- 
parison, in k' Ceti, a dwarf star with a spectra l type G5V and a 
mean rotation period of » 9 d, 'Walker et al. I (l2007 h estimated 
a relative differential rotation amplitude of ^ 0.090 + 0.006 
by modelling the optical photometry obtained by the satellite 




Time (BJD-2454900) 

Fig. 7. Amplitude of the Morlet wavelet of the total spotted area 
variation vs. the period and the time. The amplitude was normal- 
ized to its maximum value. Different colours indicate different 
relative amplitudes from the maximum (yellow) to the minimum 
(dark blue) as indicated in the colour scale in the right lower cor- 
ner. 



MOST. A result similar to that of Kepler-17 was found by 
ICroU et al.1 (l2006l) for e Eridani, a K2V star with a rotation pe- 
riod of 1 1 .4 d, showing a relative amplitude of the differential 
rotation of 0.1 1 ± 0.03. On the other hand, the more rapidly ro- 
tating star CoRoT-2 (Pmx. ~ 4.52 d) shows very little surface 
differenti al rotation with a r elative amplitude smaller than ^ 1 
percent (Lanza et al."2009a' ). In t he yo ung (^ 50 Myr) K2V 
star KIC 8429280, Frasca etaTI (l20ll found AQ/Q * 0.05 
while for the G2V stars KIC 7985 370 and KIC 7765135 , which 
have an age of 100 - 200 Myr, iFrohhch et aU (1201 2l) found 
AO/Q » 0.07 - 0.08. The differences can be attributed to the 
dependence of the differential rotation amplitude on the stellar 
rotation period and effective temperature as well as to the lim- 
ited range of latitudes covered by starspots that can vary along 
the activity cycle. Moreover, the amplitude of the differential 
rotation could change along the activity cycle in rapidly rotat- 
ing, highly active stars, such as the above mentioned Kep ler tar- 
gets that have P^pt ranging from 1.2 and 2.8 d (cf. Frasca et al. I 
l201lHFrohlichetal.ll2012h . Therefore, individual results com- 
ing from spot modelling can be compared only in a statisti- 
cal sense. Considering the general trend found bv lBarnes et al. I 
(l2005h . the amplitude of the differential rotation estimated for 
Kepler-17 appears within the expected range for a star with its 
rotation period and effective temperature. Nevertheless, a more 
extended series of data is required to derive a firm conclusion on 
this point as well as on a possi ble influence of the p lanet on spot 
activity as conjectured by, e.g.,|Lanza 1 120081 120ri1) . 

An interesting result of our analysis is the short-term spot cy- 
cle with a period of 47.1 ± 4.5 d, which is clearly detected in the 
second half of the data series. This phenomenon is reminiscent 
of the solar Rieger cycles because of its timescale and transient 
nature. The periodicity of the spotted area variations is close to 
4 times the mean rotation period of the star, while in the Sun it 
is a: 6 times and in CoRoT-2 is ss 6.5 times the rotation period. 
Those oscillations of the spotted area may be associated with 
hydromagnetic Rossby-type waves propagating in the upper part 
of the convection zone or at the interface between the radiative 
interior and the c onvection zon e where th e dynamo is proba- 
bly located (iLou I I2OOO: Zaqarashvilietal. 2010). Since only a 
few examples of stars displaying Rie ger-type cycles are known 
(cf., iMassi et al7ll2005l: iLanza et al]l2009a) . the new results on 
Kepler-17 are particularly interesting for a better understanding 
of this phenomenon in the framework of the solar-stellar connec- 
tion. 

As new Kepler data become available, it will be possible 
to refine the conclusions of the present study by investigating 
a longer time span. This is needed to assess the duration and the 
frequency of the Rieger-type oscillations in the total starspot area 
as well as to derive a definite conclusion on the possible impact 
of the planet on stellar activity. 
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